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Abstract 


When using all MSC/NASTRAN eigensolution methods except Lanczos, the analyst can replace the coupled system 
rigid-body modes calculated within DMAP module READ with mass orthogonalized and normalized rigid-body 
modes generated from the system stiffness. This option is invoked by defining MSC/NASTRAN r-set degrees -of- 
fieedom via the SUPORT Bulk Data card. The newly calculated modes are required if the rigid-body modes 
calculated by the eigensolver are not “clean” due to numerical roundoffs in the solution. When performing transient 
structural dynamic load analyses, the numerical roundoffs can result in inaccurate rigid-body accelerations which 
affect steady-state responses. Unfortunately, when using the Lanczos method and defining r-set degrees-of-freedom, 
the rigid-body modes calculated within DMAP module REIGL are retained. To overcome this limitation and to 
allow MSC/NASTRAN to handle SUPORT degrees-of-freedom identically for all eigensolvers, a DMAP Alter has 
been written which replaces Lanczos-calculated rigid-body modes with stiffness-generated rigid-body modes. The 
newly generated rigid-body modes are normalized with respect to Ihe system mass and orthogonalized using the 
Gram-Schmidt technique. This algorithm has been implemented as an enhancement to an existing coupled loads 
methodology. 
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Nomenclature 


Abbreviations 


DOF 

DMAP 

LeRC 

NASA 


F 
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X 

Y 
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a 
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a 

r 

x 
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el 

rb 


Degrees-of-freedom 

Direct Matrix Abstraction Program 

Lewis Research Center 

National Aeronautics and Space Administration 

Matrices 


Constraint forces 
Identity 
Stiffness 
Mass 

Orthogonality check matrix as computed using new r-set rigid-body modes 

Orthogonality check matrix as computed using new x-set rigid-body modes 

Cross-orthogonality check matrix 

Scaling factor for normalizing mode shapes 

Mode shape vector and matrix 

System eigenvalue 


Set Notation 


a-set (analysis DOF) 
r-set (DOF defined by SUPORT card) 
x-set (DOF in a-set having non-zero masses) 
y-set (complement of r-set in x-set) 

Superscripts 


elastic 

rigid-body 
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Introduction 


In general, to perform a modal transient loads analysis on a space structure, the system frequencies and mode shapes 
arc needed; and hence, the solution of the eigenvalue problem is required. For undamped systems, the computed 
eigenvalues correspond to the square of the natural circular frequencies (a* 2 ), and the computed eigenvector 
correspond to the normal modes of vibration (mode shapes 4>). The Givens, Householder, Inverse Power, and 
Lanczos methods are real eigensolution methods provided by MSC/NASTRAN [1]. While the Givens, Householder, 
and Inverse Power methods are executed within the DMAP module READ, the Lanczos method is executed within 
the DMAP module REIGL. Although the Givens and Householder methods are suitable for small size problems, 
they are not efficient for medium to large size problems which do not fit into the computer memory. In these cases 
a reduction procedure, such as Guyan reduction [2] or generalized dynamic reduction (GDR) [1], can be used prior 
to solving the eigenvalue problem. When combined with these reduction methods, the Givens or Householder 
eigensolvers are very efficient in solving the eigenproblem of large models. Unfortunately, they can produce poor 
results and can miss modes [3]. The Inverse Power methods are very suitable for large problems if only a small 
number of modes are required; however, they are not efficient if more than a few modes are needed [3]. The 
Lanczos method is the most accurate and efficient method that can be used to solve the eigenvalue problem for 
medium to large size models. For most aerospace applications, the finite element model of a substructure (e.g. a 
payload or a launch vehicle) is composed of a large number of DOF. Thus, the Lanczos method is the preferred 
eigensolver for these types of structures. 

Generally, aerospace systems are unconstrained (e.g. a launch vehicle flying after liftoff) and are considered free-free 
systems. For a free-free system, there are r number of rigid-body mode shapes which have zero value frequencies. 
The number of rigid-body mode shapes, r, is equal to the number of statically determinate constraints required to 
prevent all rigid-body motion. Due to numerical roundoffs, the frequencies computed by any eigensolver are not 
perfect 0.0 Hz frequencies, and the computed mode shapes are not “dean” rigid-body modes. This can result in 
inaccurate rigid-body modal accelerations which affect the steady-state responses in a transient solution. 

In order to overcome this problem, an analyst can use SUPORT Bulk Data cards with the Givens, Householder, and 
Inverse Power methods. The SUPORT card entry forces the solver to replace the computed rigid-body frequendes 
with “perfect” 0.0 Hz values. Also, the r-set DOF for the residual structure, defined via the SUPORT card, are used 
to generate rigid-body modes from the system stiffness matrix. The stiffness-generated rigid-body modes are 
orthogonalized and normalized with respect to the system mass matrix, and since they are “cleaner”, they replace 
the rigid-body modes computed from the system eigenvalue problem. Unfortunately, when using SUPORT cards 
with the Lanczos method, the rigid-body modes calculated within DMAP module REIGL are retained. The only 
effect of the SUPORT card in the Lanczos method is that it attempts to output 0.0 Hz fiequendes for the computed 
rigid-body modes without changing the originally computed eigenvectors [3]. When a SUPORT card is defined, the 
I an<y™ solver inspects the values of the computed frequendes for the first r modes. If any frequency value is near 
0.0, it is replaced by a “perfect” 0.0 value; otherwise, the computed frequency is retained. 

The objective of this work is to allow MSC/NASTRAN to handle SUPORT DOF similarly for all eigensolvers. In 
order to meet this objective, a DMAP Alter has been written which replaces Lanczos -calculated rigid-body modes 
with stiffness-generated rigid-body modes. The newly generated rigid-body modes are normalized with respect to 
the system mass and orthogonalized using the Gram-Schmidt technique. Also, the frequendes of the first r modes 
output by the Lanczos solver are replaced with “perfect” 0.0 Hz values. This DMAP Alter has been implemented 
as an enhancement to an existing coupled loads methodology at the NASA Lewis Research Center. 

The theory behind calculating system level rigid-body modes using the system stiffness, normalizing the modes with 
respect to the system mass, and orthogonal izing the modes via the Gram-Schmidt method is presented in the 
following section. In the subsequent section, implementation of the method within a DMAP Alter is described, as 
are the internal checking procedures available in the new Alter. Lastly, the new Alter is tested through a real-world 
engineering problem. The results of this example problem and the checking procedure results are presented in the 
last section. 
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Theory 


The Alter to generate system rigid-body modes using the system stiffness and enable the consistent utilization of 
SUPORT DOF with the MSC/NASTRAN Lanczos eigensolver is relatively straight forward. In Solution 63 of 
MSC/NASTRAN (normal modes for superelements), the system eigensolution is performed on the x-set DOF. These 
DOF are either the system analysis (a-set) DOF if all a-set DOF have mass associated with them, or the x-set DOF 
are a subset of the a-set DOF where only those DOF with mass associated with them have been retained. Given this, 
the coupled system eigenproblem is 


[[KJ - = {OJ 


( 1 ) 


where [h^J and [KJ are the system x-set DOF mass and stiffness matrices, and and {<|>|} are the ith system 
eigenvalue and normal mode. Assuming for a particular cutoff frequency the system has r number of rigid-body 
(rb) modes and e number of elastic (el) modes (h number of total modes), the matrix of calculated system modes, 
[*J> is written as 


i*j - [{**'> - «>*'} {♦«■} {^> - {♦*•>] 


( 2 ) 


or 


[*j - [[<■] pu] (3) 

where the first r modes, [$£ d ], are the system rigid-body modes calculated via the eigensolver. As mentioned in the 
introduction, generally these modes as calculated by the eigensolver have non-zero frequencies due to small 
numerical roundoffs. In all eigensolvers except Lanczos, these rigid-body modes are replaced by “clean” rigid-body 
modes if a SUPORT card defining r-set DOF is present in the Bulk Data deck. Hence, the system modes including 
the new rigid-body modes are written as 


py = [pc'] py] 


P) 


where [$£"] arc U»c system rigid -body modes generated using the system stiffness, normalized and orthogonalized 
with respect to the system mass. The purpose of the new Alter is to replace the rigid-body modes calculated by 
Lanczos with “clean” rigid-body modes generated using the system stiffness and corresponding to r number of r-set 
DOF defined via a SUPORT card. 

To generate new rigid-body modes with respect to the r-set DOF using the system stiffness, the system stiffness 
matrix [KJ is first partitioned according to r-set and complement (hereafter, referred to as y-set) DOF 


W [KJ 

[py py 


( 5 ) 


Next, unit motion is applied to the partitioned stiffness matrix at the r-set DOF. Since the r-set is statically 
determinate, the system should move in a rigid-body manner when unit displacements are applied sequentially to each 
r-set DOF while holding the rest and allowing the y-set DOF to move freely [4]. The system equilibrium equation 
is 
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( 6 ) 


[KJ 

py 


pj ' 


IFJ' 

py 

py 




py 


where [FJ are constraint forces required to move the system in a rigid-body manner (these forces should be 
computed zeros). From the lower partition of Eq. ( 6 ), the un-nonnalized, un-orthogonalized rigid-body modes for 
the y-set DOF are calculated as constraint modes [4] 


[*7i - -py'py 


(7) 


Hence, the un-normalized and un-orthogonalized x-set system rigid-body modes are 


PCI 


[!„] 

[*7i 


( 8 ) 


where [y is an identity matrix of order r which represents the unit translations (and/or rotations) applied at the r-set 
DOF to generate the rigid-body modes. 


To orthogonalize and normalize the computed rigid-body modes with respect to the mass matrix, first the system 
rigid-body mass matrix is calculated as 


[m„ a ] = [$rn M jr*;n 


(9) 


For each rigid-body mode, the approach taken is to first orthogonalize the mode to all previous rigid-body modes 
and then normalize it with respect to the system mass given in Eq. (9). Obviously, the first rigid-body mode needs 
only to be mass normalized. All other modes are orthogonalized with respect to the previously calculated modes. 
The procedure for orthogonalizing and normalizing is performed on the r-set partition of the rigid-body modes. This 
is much more efficient than performing the operations on the foil x-set DOF modes. After orthogonalization and 
normalization of the r-set DOF modes, the rigid-body modes are expanded to the x-set DOF. 

Let the r-set partition of the un-orthogonalized and un-normalized system rigid-body mode shapes of Eq. ( 8 ) (an 
identity matrix) be written as 


i*n = pj = [{«\> «?} ~ <♦:>] 

For the first rigid-body mode, only mass normalization is required. Hence, 


{4>r> - 


where the scaling factor 04 is calculated as 


a i 


1 


( 10 ) 


(u) 


( 12 ) 
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Therefore, {<f>j} is the first rigid-body mode normalized with respect to the system mass. 

For each of the remaining r-2 number of rigid-body modes, the following procedure is performed. For a given jth 
mode, the mode is first orthogonalized to the previous j-1 number of modes using the Giam-Schmidt 
orthogonalization scheme [5], or 


{♦!> - - E({^} T [M B *]{4.! E }){^} 

k-1 


(13) 


After orthogonalization of the jth mode, it is normalized with respect to the system mass as 

m = <*.{♦!> 


where the scaling factor is calculated as 


= 1 


(14) 


(15) 


Given that all system rigid-body modes have been orthogonalized and normalized with respect to the system mass 
for all r-set DOF, the final modes corresponding to the r-set DOF are written as 

[*n = [{♦?> (♦'} - {♦;}] (16) 

This matrix can be considered a scaling matrix required to linearly transform the rigid-body modes given by Eq. (8) 
to form the final orthogonalized and mass normalized x-set DOF rigid-body modes. The final new system rigid-body 
modes expanded to the x-set DOF are 


m = [*7][*r] 


i*n 

Nn 


where 


[*;n - 


(17) 


(18) 


To prove that the final computed x-set rigid-body modes are mass orthogonalized and normalized, consider the 
following triple matrix multiplication 


[*ri T [mj m 


Substituting Eq. (17) into (19) 

[*7Y [mj [«n - (i*ri [*n) T mj i*n) 

- [<i»r] T (f^TTWJ[*n) 1®T) 


(19) 


( 20 ) 
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The middle term between parentheses in Eq. (20) is the rigid-body mass matrix given by Eq. (9), hence 


[*ri T [m j [*n - [*n T [*ci i«n 


( 21 ) 


Since the r-set rigid-body modes given by Eq. (16) were orthogonal! zed and normalized with respect to the rigid- 
body mass matrix using the Gram-Schmidt orthogonalization scheme, the right hand side of Eq. (21) is equal to an 
identity matrix. Therefore Eq. (21) becomes: 

[*rr [mj [*n - iu (22) 


This proves the final system x-set rigid-body modes given by Eq. (17) are mass orthogonalized and normalized. 


Implementation 


In order to replace system rigid-body modes calculated using the Lanczos eigensolver with rigid-body modes 
generated using the system stiffness, a new DMAP Alter has been implemented in the NASA LeRC coupled loads 
methodology. This Alter to MSC/NASTRAN Solution 63 is included in the Appendix. 

The following are the alterations made to Solution 63: 

• Prepare required data blocks and replace the solved-for rigid-body eigenvalues with perfect zeros. 
This is performed first by transforming the LAMA table to a matrix and then replacing the rigid- 
body frequency entries with zero values. The discarded rigid-body frequencies output by Lanczos 
are printed for inspection. 

• Use residual USET table to form a partition vector, VXCOMPRN, for partitioning out r-set DOF 
from x-set DOF. 

• Calculate system rigid-body modes from system stiffness as given by Eqs. (7) and (8). Generate 
system rigid-body mass matrix, as defined by Eq. (9), using new modes and print as a check. 

• Normalize first new system rigid-body mode with respect to system mass as shown in Eqs. (11) 
and (12). 

• Use the DO WHILE DMAP statement to loop over the rest of the new system rigid-body modes. 
Orthogonalize each subsequent mode with all previous modes using Gram-Schmidt method (Eq. 
(13)) and normalize each mode with respect to system mass (Eqs. (14) and (15)). 

• Expand new system rigid-body modes from r-set DOF to x-set DOF (Eqs. (17) and (18)). 

• Check orthogonality of new system rigid-body modes if required by the analyst via new parameter 
CHECKRBM. The checks performed are for the r-set DOF rigid-body modes: 

= [XJ (23) 


and for the expanded x-set DOF rigid-body modes: 


[«STTWjm - Py 


(24) 
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where [Xj and [YJ matrices are printed for inspection. Theoretically, these matrices are equal 
to the identity matrix [IJ. 

• Replace old system rigid-body modes in PHiX datablock, generated by REIGL DMAP module, 

with new system rigid-body modes. 

With that, the system modes are ready for further operations within Solution 63. 


Numerical Example 

Fora launch vehicle/spacecraft coupled system, rigid-body mode shapes were computed using three different methods 
with the same SUPORT DOF defined. The three methods used were Lanezos, Inverse Power, and Lanczos with the 
new Alter. The order of the x-sized matrices for which the eigensolutions were performed was equal to 15,616. 
The purpose of solving the eigenproblem with three different methods was to illustrate the following: 

1- Rigid-body mode shapes output by the Lanczos solver (without the new Alter) are different than 
those output by other solvers (Inverse Power without the new Alter) when the same SUPORT DOF 
are defined. 

2- Stiffness -generated mass orthogonalized and normalized rigid-body modes are accurately computed 
by the new Alter when used with the Lanczos solver with SUPORT DOF defined. 

3- Rigid-body mode shapes output by the Lanczos solver (with the new Alter) are similar to those 
output by other solvers (Inverse Power without the new Alter) when the same SUPORT DOF arc 
defined. 


To show that the rigid-body mode shapes computed by Lanczos (without the Alter) arc different than those computed 
by Inverse Power (without the Alter) with the same SUPORT DOF defined, cross-orthogonality checks were used. 
The cross-orthogonality matrix [W n ] was computed as: 

[M^] _ [^j (25) 


where the SINV superscript denotes the Inverse Power method. As shown in Table (1), the result of this check is 
not an identity matrix which demonstrates the difference between the rigid-body modes output by both eigensolvers. 

To determine whether or not the new Alter properly mass normalized and orthogonalized the SUPORTed rigid-body 
modes, the orthogonality checks shown by Eqs. (19) and (20) were performed. Results were excellent and are shown 
in Tables (2) and (3). 

To determine whether or not the Alter generated SUPORTed rigid-body modes (Lanczos with Alter) equal to those 
generated using standard MSC/NASTRAN (Inverse Power without Alter), the former were compared with the latter 
via cross-orthogonality checks and direct comparisons. The cross-orthogonality matrix was computed as: 

[•TFmj ■ [ Z J C 26 ) 


The result of this check, shown in Table (4), is an identity matrix which shows the Alter property generates 
SUPORTed rigid-body modes as compared to standard MSC/NASTRAN. The rigid-body modes were then 
subtracted from one another in a direct comparison as 

[<j)^ INV ] _ [<j)^ MCOTWitbo,rtAlte ] s [oj (27) 
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Aga in, the results were excellent, where the largest term in the difference was on the order of 10 1 (calculated zero). 
Given these comparison results on a real-world engineering problem, it has been shown that the new Alter performs 
as required and replicates the operations performed by MSC/NASTRAN for all other eigensolution methods when 
SUPORT DOF are defined. 


Conclusion 


In order to enable MSC/NASTRAN to handle SUPORT DOF with the Lanczos method in a similar manner to the 
other eigensolvers, a Solution 63 DMAP Alter has been written as an enhancement to the NASA LeRC coupled loads 
methodology. The new Alter replaces the rigid-body modes generated by DMAP module REIGL with rigid-body 
modes computed from the system stiffness matrix. The newly generated rigid-body modes are orthogonalized and 
normalized with respect to the system mass matrix using the Gram-Schmrdt orthogonalization scheme. In addition, 
the rigid-body frequencies output by the Lanczos solver are replaced with perfect 0.0 Hz values. A real-world 
engineering problem was used to test the accuracy of the newly developed Alter. Rigid-body modes were calculated 
using Lanczos (with the Alter) and Inverse Power (without the Alter) methods with the same SUPORT DOF defined. 
It has been shown, through a comparison of the generated rigid-body modes, that using the Alter with the Lanczos 
method produces very accurate results. 
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-6.9215e-01 2.0783e-03 5.8888e-03 1.9600e-03 7.2172e-01 -9.1999e-04 

5.941 le -04 6.8867e-01 -2.4478e-03 7.2361e-01 -3.4169e-03 -43826e-02 

7.6239e-03 1.5092e-03 9.9997e-01 1.86 13c -03 -8 .585 le -04 -1.1783e-03 

4.1390e-02 2.0442c-02 7.8191e-04 43854e-02 4.0781e-02 9.9714c-01 

-7.2033c -01 1.6575c -02 4.9701c-03 -1.4726c-02 -6.9079c-01 5.8456e-02 

1.6714c -02 7.2459c -01 9.1000c-05 -63865c-01 13830c -02 1.4091e-02 


Table 1. -- Cross-orthogonality check computed using SUPORTed rigid-body modes 
generated using Lanczos without the new Alter vs. Inverse Power 


py 


l.OOOOc+OO 
1.0842c-19 
-2.7105e-20 
-2.1684e-19 
-3.5527e-15 
-8.3267c- 17 


1.0842c-19 

1.0000e+00 

0.0000e+00 

0.0000e+00 

-5.4210e-20 

1.1102c-16 


-2.7105e-20 
O.OOOOc+OO 
1.0000c +00 
-1.0408c-17 
5.5511e-17 
6.1240c-19 


-2. 1684c -19 
0.0000e+00 
-1.0408c-17 
1.0000c +00 
1.3553c-19 
-2.2204c-16 


-3.5527c -15 
-5.4210e-20 
5.5511e-17 
1.3553c-19 
1.0000c +00 
5.551 le -17 


-8.3267c -17 
1.1102c-16 
6.1240e-19 
-2.2204c -16 
5.5511c-17 
1.0000c +00 


Table 2. — Orthogonality check computed using r-set size 
SUPORTed rigid-body modes generated using new Alter. 


1.0000c +00 -1.4993e-14 
-1.4993e-14 1.0000c +00 


-5.0910c-16 

3.8540e-15 

-2.1588c-12 

-1.3701c-14 


-7.1370e-17 
7.003 lc-13 
6.0921e-14 
2.8727e-14 


-5.0910e-16 
-7.1370e-17 
1.0000c +00 
2.6143C-15 
1.6721c-15 
5.4753c-15 


3.8540c-15 -2.1588c-12 
7.003 lc-13 6.0921c-14 

2.6143C-15 1.6721c-15 

l.OOOOc+OO -4.0034c -14 
-4.0034c-14 l.OOOOc+OO 
1.9734c-14 9.0601c-14 


-1.3701c-14 
2.8727c-14 
5.4753c-15 
1.973 4c-14 
9.060 lc-14 
l.OOOOc+OO 


Table 3. — Orthogonality check computed using x-set size 
SUPORTed rigid-body modes generated using new Alter. 
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1.0000e+00 

3.0538e-15 

-2.7188e-16 

-6.2457e-15 

-3.9616e-13 

3.6614c-14 


1.1332e-14 

l.OOOOc+OO 

-1.6059e-16 

3.9230e-14 

-7.0599e-14 

-3.6741e-14 


-9.5581e-17 
-5.1448e-16 
1.0000e+00 
-2.4279c -16 
-2.5558c -15 
2.2042e-16 


-1.6161c-14 
75368c-13 
1.0776e-15 
1.0000c +00 
1.8957e-14 
-8.1182c-15 


-2.0192c -12 
1.1361e-13 
-2.4310e-16 
-1.0624c-13 
1.0000c +00 
4.4146c-14 


-2.1386e-14 
-1.4128c-14 
53 789c-15 
3.9857c-14 
3.7061e-14 
1.0000c +00 


Table 4. -- Cross -orthogonality check computed using SUPORTed rigid-body modes 
generated using Lanczos with the new Alter vs. Inverse Power 


Appendix 


$ LANCZOS E I GEN SOLUTI ON WITH PROPER SUPORT OF R-B MODES (RF63D401) 

$ 

$ RIGID FORMAT 63 - NORMAL MODES WITH SUPERELEMENTS 

$ MSC/NASTRAN VERSION 67 

$ 

$ 

$ THIS DMAP ALTER IS USED TO GENERATE RIGID-BODY MODES USING SUPORT 
$ DOF WITH THE LANCZOS EIGENSOLUTION METHOD. THE CALCULATED RIGID- 
$ BODY MODES GENERATED WITHIN REIGL ARE DISCARDED AND REPLACED WITH 
S RIGID-BODY MODES GENERATED FROM THE SYSTEM STIFFNESS. THESE NEW 
$ MODES ARE MASS NORMALIZED AND ORTHOGONALIZED USING THE GRAM-SCHMIDT 
$ METHOD . 

$ 

$ 

$ REQUIREMENTS TO USE THIS DMAP SEQUENCE - 

$ 

$ FILE MANAGEMENT SECTION - 

s 

$ NO SPECIAL REQUIREMENTS 

$ 

$ EXECUTIVE CONTROL DECK - 

$ 

$ DIAG 8,14,20 RECOMMENDED 

$ SOL 63 

$ COMPILE SOL63,SOUIN=MSCSOU, NOLIST, NOREF 

$ INCLUDE THIS ALTER 

$ 

§ CASE CONTROL DECK - 

$ 

$ NO SPECIAL REQUIREMENTS 

$ 

$ BULK DATA DECK - 

$ 

$ THE FOLLOWING PARAMETERS MAY BE DEFINED - 

$ 

$ PARAM, CHECKRBM - IF .GE. 0, RIGID-BODY MASS AND ORTHOGONALITY 

$ (OPTIONAL) MATRICES ARE GENERATED AND PRINTED AND CAN BE USED 

$ TO CHECK THE NEW RIGID-BODY MODES. 

$ (DEFAULT = -1 - DO NOT PERFORM CHECKS) 

$ 

$ EXAMPLE NASTRAN DECK - 

$ 

$ ID LANCZOS, SUPORT 

$ SOL 63 

$ TIME 30 

$ DIAG 8,14,20 

$ COMPILE SOL63,SOUIN=MSCSOU, NOLIST, NOREF 

$ INCLUDE RF63D401 . V67 

$ CEND 

$ TITLE = SYSTEM MODES 

$ SUBTITLE = LANCZOS WITH SUPORT CARDS 

$ $ 

$ SET 1000 = 0 

$ SEALL = 1000 

$ $ 

$ DISP(PLOT) = ALL 

§ $ 

$ SUBCASE 1 

$ SUPER = 200 

$ LABEL = SPACECRAFT 

$ METHOD = 100 
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$ SUBCASE 2 

$ SUPER - 100 

$ LABEL = LAUNCH VEHICLE 

$ METHOD = 100 

$ SUBCASE 3 

$ $ SUPER = 0 

$ LABEL = COUPLED SYSTEM 

$ METHOD =50 

§ $ 

$ BEGIN BULK 

s 

$ . STANDARD BULK DATA INCLUDING SUPORT CARDS FOR RIGID-BODY MODES 

$ 

s $ 

$ $ perform checks on new rigid-body modes 

$ § 

$ PARAM, CHECKRBM, +1 

$ $ 

$ ENDDATA 

$ 

$ HISTORY DOCUMENTATION - 

$ 

$ VERSION 1.0 01 -OCT- 9 3 A. ABDALLAH, A. BARNETT, T. WIDRICK 

$ - ORIGINAL VERSION 

$======================*«==«=========================================== 

$23456789012345678901234567890123456789012345678901234567890123456789012 

$1234567 

$ 

ALTER 982,982 $ V67 REPLACE REIGL 

$ 

$ DEFINE DEFAULT VALUES 
$ 

TYPE PARM, , I , Y, CHECKRBM=-1 $ DO NOT PRINT CHECKS 

$ 

$ RENAME LAMA AND PHIX FROM REIGL 

REIGL KXX, MXX, DYNAMICS, CASES, ,MR,DMX, VXCOMPR/ 

LAMAO,PHIXO,MI,EIGVMAT,OUTVEC/V,N,READAPP/S,N,NEIGV $ 


GENERATE RIGID-BODY AND ELASTIC MODES PARTITIONING VECTOR 


$ 

$ 

$ 

LAMX 

MATMOD 

MATMOD 

PARAML 

MATGEN 

$ 

$ SET 
$ 

PARTN 
MESSAGE 
MESSAGE 
MESSAGE 
MATPRN 
MERGE 
MATGEN 
MATGEN 
MERGE 
MERGE 
TRNSP 
LAMX 
§ 

ALTER 
$ 


, , LAMAO/LAMAOM/ - 1 $ 

LAMAOM, , , , , /FREQ, /1/3 $ 
LAMAOM, , , , , / GNMASS , / 1 /4 $ 
LAMAOM// ' TRAILER ' /2/S,N, NFREQ 
, /VHRBEL/6/NFREQ/NORSET/NFREQ 

RIGID-BODY FREQUENCIES TO ZERO 


FREQ, ,VHRBEL/FREQRB,FREQEL, , /I $ 

//' '/ $ 

//'RIGID-BODY FREQUENCIES AS OUTPUT BY LANCZOS'/ $ 

//' '/ $ 

FREQRB// $ 

, , FREQEL , , , , VHRBEL /FREQEL2 / 1 $ 

, / VLAM1 /4/1/2//1//2 $ 

, / VLAM2 /4/1/3//1//3 $ 

FREQEL2 , , , , VLAM1 , /LAMA1 / 1 $ 

LAMA1,, GNMASS,, VLAM2 , /LAMA2/1 $ 

LAMA2/LAMA3 $ 

LAMA3, /LAMA $ 


FORM MATRIX FROM LAMAO 
EXTRACT FREQUENCIES 
EXTRACT GEN MASSES 
NFREQ = NO. OF FREQ 
H = RB / EL 


PARTITION RB/EL FREQ 


PRINT LANCZOS RB FREQ 
ROW MERGE WITH NULLS 
1 

0 1 

WITH NULLS 
W/ GEN MASS 


VLAM1 = 0 
VLAM2 = 0 
COL MERGE 
COL MERGE 
TRANSPOSE 
FORM NEW LAMA 


TABLE 


984 $ V67 


AFTER OFP OF LAMA 
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REPLACE RIGID-BODY MODES IN PHIX BY RIGID-BODY MODES GENERATED 
FROM STIFFNESS MATRIX ABOUT THE R-SET DOF. THE RIGID-BODY MODES 
ARE NORMALIZED AND ORTHOGONALIZED WITH RESPECT TO THE MASS MATRIX 
USING THE GRAM-SCHMIDT ORTHOGONALIZATION SCHEME. 


V6 = OVRWRT $ 


$ 

$ 
s 
$ 

$ 

FILE 

$ 

$ GENERATE R-SET AND COMPLEMENT OF R-SET PARTITIONING VECTOR 

$ 

IF (NOARED = -1) THEN $ 

VEC USET/VXCOMPRN/'A'/'COMP'/'R' $ 

ELSE $ 

PARTN VACOMPR, , VAXW/VXCOMPRN, , ,/l $ 

ENDIF $ 

$ 

$ CALCULATE RIGID-BODY MODES FROM KXX 

$ , * 

PARTN KXX, VXCOMPRN, /KCRCR, ,KCRRF , /— 1 $ 

SOLVE KCRCR,KCRRF/DMCR//-1 $ 

MATGEN , /DMRF/l/NORSET $ 

MERGE DMCR , DMRF , , , , VXCOMPRN /DMXN / 1 $ 

SMPYAD DMXN, MXX, DMXN, , , /MRIG/3////1////6 $ 

MESSAGE //''/$ 

MESSAGE //'SYSTEM RIGID-BODY MASS MATRIX ABOUT R-SET DOF'/ $ 
MESSAGE //''/$ 

MATPRT MRIG/ / $ PRINT MRIG 

$ 

$ NORMALIZE FIRST RIGID-BODY MODE (OPERATING ON R-SET ONLY) 

$ 

MATGEN , /VFIRST/6/NORSET/1/NORSET $ 

PARTN DMRF,VFIRST, /INPHI1 , , , / I $ 

SMPYAD INPHI 1 , MRIG, INPHI 1 , , , /GNM1/3////1////6 

PARAML GNM1//'DMI' /1/1/S, N,GNM2 $ 

PARAMR // ' SQRT ' /S , N , GNM3 /GNM2 $ 

PARAMR / / ' DIV ' /S , N , GNM4 / 1 . 0 / GNM3 $ 

PARAMR / / ' COMPLEX ' / /GNM4 /0.0/S,N, GNM4CX S 

ADD INPHI 1 , /INPHI2/GNM4CX $ 

MERGE INPHI2, ,, ,VFIRST, /PHI/1 $ 

PARAM / / ' SUB ' /S,N, VCOUNT/NORSET/2 $ 

MATGEN , /V6/7/NORSET/NORSET/1 $ 


ALLOW OVERWRITE 


A - COMP / R 
A = COMP / R 


SYMMETRIC PARTITION 
FORM RIGID-BODY TRANS 
GENERATE IDENTITY MTX 
ROW MERGE 

MRIG = RIGID-BODY MASS 


ZERO FOR 1ST MODE 
COLUMN PARTITION 
$ GNM1= GEN MASS 1ST MD 
CHANGE TO SCALAR 
GNM3 = SQRT ( GNM2 ) 

GNM4 = INV ( GNM3 ) 

CHANGE TO COMPLEX 
MASS NORMALIZE 
COLUMN MERGE W/ NULLS 
VCOUNT = NORSET - 2 
GENERATE NULL MATRIX 


$ 

$ MASS NORMALIZE AND ORTHOGONALIZE EACH SUCCESSIVE RIGID-BODY MODE 
$ WITH THE PREVIOUS MODES USING GRAM-SCHMIDT METHOD 
$ 

(VCOUNT >= 0) $ 

//' SUB '/S,N,INDX/NORSET /VCOUNT $ 

/ / ' SUB ' /S , N , INDX1 / INDX/ 1 $ 

, /VINDX/6/NORSET/0/INDX1/1/VCOUNT $ 

DMRF , VINDX ,/Vl,,,/l $ 

PHI, MRIG, VI, , , /V2/3/ ///l $ 

PHI,V2,V1/V3//-1 $ 

V3 , MRIG , V3 , , , / SMI /3////1////6 $ 

SMI// 'DMI' /1/1/S, N,SM2 $ 

//'SQRT'/S,N,SM3/SM2 $ 

/ / ' DIV ' /S , N , SM4 / 1 . 0 / SM3 $ 

// 'COMPLEX' //SM4/0. 0/S, N,SM4CX $ 

V3 , /V4 / SM4CX $ 

V4, , , , VINDX, /V5/1 $ 

PHI /V6 /NEVER $ 

PHI , V5 /V6 $ 

V6 /PHI /ALWAYS $ 

= VCOUNT - 1 $ 


DO WHILE 
PARAM 
PARAM 
MATGEN 
PARTN 
SMPYAD 
MPYAD 
SMPYAD 
PARAML 
PARAMR 
PARAMR 
PARAMR 
ADD 
MERGE 
EQUIVX 
ADD 
EQUIVX 
VCOUNT 

ENDDO $ 

S 


TOP OF LOOP 

INDX = NORSET - VCOUNT 

INDX1 = INDX - 1 

GEN PARTITION VECTOR 

COLUMN PARTITION 

V2 = (PHI)T * MRIG * VI 

V3 = VI - PHI*V2 

SMI - (V3)T * MRIG * V3 

CHANGE TO SCALAR 

SM3 - SQRT ( SM2 ) 

SM4 = INV ( SM3 ) 

CHANGE TO COMPLEX 
MASS NORMALIZE 
COL MERGE WITH NULLS 
BREAK EQUIVALENCE 
V6 = PHI + V5 
MAKE EQUIVALENT 
RESET VCOUNT 
BOTTOM OF LOOP 
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$ EXPAND RIGID-BODY MODES FROM R-SET TO X-SET 

$ 

MPYAD DMCR, PHI , /PHICR $ PHICR = DMCR * PHI 

MERGE PHICR, PHI, ,, , VXCOMPRN/DMRGFIN/1 $ ROW MERGE 

S 

$ CHECK ORTHOGONALITY OF COMPUTED RIGID-BODY MODES IF REQUIRED BY 
$ PARAM , CHECKRBM 
$ 

IF (CHECKRBM > -1) THEN $ 

MESSAGE //' '/ $ 

MESSAGE //'RIGID-BODY SCALING MATRIX TO ORTHOGONALIZE R-B MODES'/ $ 
MESSAGE //' '/ $ 

MATPRT PHI// $ PRINT RB MD SCALING MTX 

SMPYAD PHI ,MRIG, PHI, , , /MRBRSET/3// / /l/ ///6 $ MRBRSET= ( PHI ) T*MRIG*PHI 

MESSAGE //' '/ $ 

MESSAGE //'ORTHOGONALIZED RIGID-BODY MASS MATRIX COMPUTED USING'/ $ 
MESSAGE //' R-SET RIGID-BODY MODES AND MASS'/ $ 

MESSAGE //' '/ $ 

MATPRT MRBRSET// $ PRINT RB MASS MATRIX 

SMPYAD DMRGFIN,MXX,DMRGFIN, , , /MRBXSET/3////1////6 $ TRIPLE MULTIPLY 
MESSAGE //' '/ $ 

MESSAGE //'ORTHOGONALIZED RIGID-BODY MASS MATRIX COMPUTED USING'/ $ 
MESSAGE //' X-SET RIGID-BODY MODES AND MASS'/ $ 

MESSAGE //' '/ $ 

MATPRT MRBXSET// $ PRINT RB MASS MATRIX 

ENDIF $ 

S 

$ REPLACE RIGID-BODY MODES IN PHIXO WITH COMPUTED RIGID-BODY MODES 

$ 

PARTN PHIXO, VHRBEL, /, , PHIXOEL, /I $ COLUMN PARTITION 

MERGE DMRGFIN, , PHIXOEL, , VHRBEL, /PHIX/1 $ COLUMN MERGE 
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